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We study a cellular automaton opinion formation model of Ising type, 
with antiferromagnetic pair interactions modeling anticonformism, and fer- 
romagnetic plaquette terms modeling the social norm constraints. For a 
sufficiently large connectivity, the mean-field equation for the average mag- 
netization (opinion density) is chaotic. This "chaoticity" would imply ir- 
regular coherent oscillations of the whole society, that may eventually lead 
to a sudden jump into an absorbing state, if present. 

However, simulations on regular one-dimensional lattices show a dif- 
ferent scenario: local patches may oscillate following the mean-field de- 
scription, but these oscillations arc not correlated spatially, so the average 
magnetization fluctuates around zero (average opinion near one half). The 
system is chaotic, but in a microscopic sense where local fluctuations tend 
to compensate each other. 

By varying the long-range rewiring of links, we trigger a small-world 
effect. We observe a bifurcation diagram for the magnetization, with period 
doubling cascades ending in a chaotic phase. As far as we know, this is the 
first observation of a small-world induced bifurcation diagram. 
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The social implications of this transition are also interesting. In the 
presence of strong "anticonformistic" (or "antinorm" ) behavior, efforts for 
promoting social homogenization may trigger violent oscillations. 

PACS numbers: 05.45.Ac,05.50.+q,64.60.aq,64.60.Ht 

1. Introduction 

Social norms are the basis of a community. Social norms are often 
adopted and respected even if in contrast with an individual's immediate 
advantage, or, alternatively, even if they are costly with respect to a "naive" 
behavior. Indeed, the social pressure towards a widespread social norm is 
sometimes more powerful than a norm imposed by punishments. 

On the other hand, it is well known that the establishment of social 
norms is difficult to plan, and their imposition is hard to be fulfilled. This 
problem has been affronted by Axelrod in a game-theoretic formulation [I], 
as the foundation of the cooperation and of the society itself. Axelrod's 
idea is that of a repeated game. Although in an one-shot game it is always 
profitable to win not following any norm, in a repeated game there might 
be several reasons for cooperation j2], the most common ones are direct 
reciprocity and reputation. In all these games, the crucial parameters are 
the cost for cooperation with respect to defeat, and the expected number 
of re-encounters with one's opponent or the probability that one's behavior 
will become public. One can assume that these aspects are related to the 
size of the local community with which one interacts, and on the fraction 
of people in this community that share the acceptance of the social norm. 
Indeed, the behavior of a spatial social game is strongly influenced by the 
network structure [3] . 

The previous approach assumes perfect rationality of agents, and does 
not take into consideration "irrational" tendencies like for instance educa- 
tion. It is well known that a given predisposition towards conformism or 
anticonformism (i.e., the education by parents, school and the social com- 
munity) may influence the acceptance of a given social norm. 

In this paper we model the dynamics of the acceptance of a social norm 
in a community with different degrees of conformism or anticonformism. 
We shall consider a simplified cellular automata model, already introduced 
in Ref. jl]. See Ref. [5] for a detailed review of models of social dynamics. 

2. Cellular automata model 

Let us denote by s« = Si(t) the opinion of individual i at time t. We 
consider the case of two opinions in competition Si £ {0, 1}. We can switch 
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to Ising-like variables (spin) Oi G {—1,1} by the transformation 

Oi = 2si - 1. 

The individual opinion evolves in time according to the opinions of neigh- 
bors, identified by an adjacency matrix a^j £ {0, 1}. This matrix defines 
the network of interactions and is considered fixed in time. An individual 
may be part of his neighborhood. 

The neighborhood of an individual i is the set of individuals j such that 
dj j = 1 . The connectivity ki of individual i is the number of nonzero entries 
in ctjj, i.e., 

ki = a i,j- 
3 

In the following, we shall consider uniform neighborhoods with connec- 
tivity k, where the actual neighborhood can be either regular 



1 if j — i < k, 
otherwise. 



or partially rewired, where a fraction p of the nonzero entries of each row of 
cijj is set to zero, and replaced with a randomly chosen j' (with djj/ = 0), 
setting djj/ = 1. 

We choose to assign equal weight to all neighbors, so we define the local 
field (social pressure) hi as 

hi = . (1) 

The local field takes values between and 1. We might also add an external 
field H, modeling written or broadcasting media, but in this study we always 
keep H = 0. 

We are modelling here a completely uniform society, i.e., we assume 
that the individual variations in the response to stumuli are quite small. 
Moreover, we do not include any memory effect, so that the dynamics is 
completely Markovian. 

The effects of the social pressure is assumed to be proportional to the 
social field. The evolution is defined by the transition probabilities 

r(l\hi), 

denoting the probability of observing a spin S{(t + 1) = 1 given a local field 
hi at time t. Clearly, r(0|/ij) = 1 — r(l\hi). 
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By denoting by s a spin configuration (si, s 2 , • • • ), and by P(s, t) the 
probability of observing it at time t, we have 

P(s',t + l) = Y,W(s'\s)P(s,t) (2) 

s 

with 

w( S 'i S )=n r (*)' 

i 

and hi given by Eq. ([I]) . 

If no transition probability is zero or one, we can map the Markov matrix 
W onto a dynamic equilibrium model, of Ising type [6] . 

Eq. (|2j may be equivalently written using restricted distributions. Let us 
call p n (si,..., Si +n ; t) the (restricted) probability of observing the sequence 
Si, ... , Si +n at time t, defined as 

Pn(Si > • • • ? Si+n j t) — ^ ^ P^Sl j • • • > Sj , . . . , Sj_|_ n , . . . , t) . 

s J :j^{i,....i+ra} 

By assuming spatial uniformity, the p n do not depend on the index i. The 
quantity pi(l,i) is the usual density, that will be denoted by c. 

Let us consider as an example the one-dimensional regular lattice case 
with uniform connectivity k = 2, an = a^j+i = 1. Using the restricted 
distributions, the evolution equation of the system is given by the infinite 
hierarchy 

pi(s' 1 ;f + l) = ^2 P2(si,s 2 ;t)T(s' 1 \s 1 ,s 2 ) 
P2(s' 1 ,s' 2 ;t + 1) = ^2 P3{si,S2,s 3 ;t)T(s' 1 \si,s 2 )T(s'2\s 2 ,s 3 ) (3) 

Sl,S2iS3 

p 3 (si,s2)4;* + 1 ) 



where for readability (and generality) we have written t^s'^si, s 2 ) instead 
afr(si|(fli + « 2 )/2). 

3. Phase transitions in uniform societies 

Let us express the transition probabilities as 

if h < q, 



r(l\h) 



' if q < h < 1 - q, (4) 



1 + exp(-2 J(2h - 1)) 
1 - e if h > 1 - q 
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shown in Fig. [T] In this way we model a standard dynamic Ising model with 
ferromagnetic (for J > 0) or antiferromagnetic (J < 0) interactions, with 
plaquette terms given by £q For e = (infinite plaquette terms) we have 
the absorbing states s = and s = 1 if the social pressure is above (below) 
the threshold 1 — q (q), respectively. 

In one dimension, with k = 3, 1/3 < q < 1/2 and e = this model 
exhibits a nontrivial phase diagram [7] , with two directed-percolation tran- 
sition lines that meet a first-order transition line in a tricritical point, be- 
longing to the parity conservation universality class. Essentially, we have 
the stability of the two absorbing states for J > (ferromagnetic interac- 
tions, conformistic society, ordered phase), while for J < (anti-ferro and 
anti-conformistic) the absorbing states are unstable and a new, disordered 
active phase is observed. This scenario corresponds essentially to the sim- 
plest mean-field picture. The interpretation of mean-field predictions is that 
< c < 1 corresponds to the active phase, which is microscopically chaotic, 
with the appearance of transient correlated patches ("triangles"). This is 
due to the presence of the unstable absorbing states: occasionally a patch of 
sites "fall" into one of these states. Since the absorbing state is absorbing, 
it can be abandoned only by "erosion" at boundaries, and this originates 
the "triangular" pattern. 

The model has been studied also in the one- dimensional case with larger 
neighborhood [8]. In this case we observe again the transition from an or- 
dered to an active, microscopically chaotic phase, but this transition occurs 
through a disordered phase, with no apparent structure in the time-space 
pattern, which is moreover insensitive to variations of J. Indeed, if the 
system "falls" into a truly disordered configuration, then h is everywhere 
equal to 0.5 and the transition probabilities r become insensitive to J and 
equal to 0.5, Eq. Q. This disordered regime is therefore different from the 
"microscopically chaotic" one. 



4. Mean-field chaotic behavior 

The mean field approximation consists in truncating the previous hier- 
archy, Eq. (|3]), at a given point by factorizing the distribution probabilities. 
At the lowest level, for a uniform connectivity k, one gets 

c = /(c) = £ ( k ) c*(l - c) k - v r{l\v/k\ (5) 

which is also shown in Fig. [Tj In this approximation, the return map Eq. (|5]) 
may become chaotic, for antiferromagnetic (J < 0) couplings and sufficiently 



1 A similar model can be defined using a standard Hamiltonian formalism by including 
ferromagnetic plaquette terms proportional to (2/i — l) 3 or higher odd powers. 
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9 0.2 0.4 0.6 0.8 1 -9 1 

h, c 

Fig. 1. Transition probabilities r(l\h) (Eq. Q, thick dashed line) and return map 
/(c) (Eq. (J5j, thick continuous line) as a function of c. Some iterations of the map 
are also shown (thin continuous line) and the bisectrix (dash-dotted line). The 
threshold values q and 1 — q are marked by thin dashed vertical lines, and the 
values e and 1 — e by thin dashed horizontal lines. Here k — 20, q = 0.1, J = —3, 
e = 0.2. 




1 1 1 1 1 1 1 

1 2 3 4 5 6 

-J 

Fig. 2. Mean field bifurcation diagram of Eq. (|) for < —J < 6 and k = 20, 
q = 0.1, e = 0.2. Transient of 100 time steps, and plot of 10 iterations for 4 random 
initial conditions. 



large neighbors, see Fig.[T] Notice that in this case, we do not have absorbing 
states, but the plaquette terms are necessary to give origin to the chaotic 
oscillations. 

By varying J, one can observe a bifurcation diagram of "logistic" type 
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Fig. 3. Mean field bifurcation diagram ofEq. ^ for J = -6 and varying k, q = 0.1, 
e = 0.2. Transient of 100 time steps, and plot of 10 iterations for 20 random initial 
conditions. 



(period doubling) as shown in Fig. [2j 

The connectivity k plays a fundamental role. As shown in Fig. [3] for 
J = —6, chaotic oscillations occur only for some values of k. For J = —6, 
k = 20 corresponds to an almost compact chaotic region, near to a window 
with a strongly periodic behavior. 

Since we have approximated the behavior of an extended system with a 
scalar equation, we have imposed spatial homogeneity, which is not generally 
observed. Better approximations are obtained by replacing Eq. §5§ with a 
spatial, coarse-grained description as follows 

c(t + l) = f(c) + D^+ V ^rrcj, (6) 

where x represent the coarse-grained space index, c = c(x,t), D a diffusion 
coefficient, that also plays the role of a surface tension term that tends to 
make the system homogeneous. The term r/^ c(l — c) represents the local 
fluctuations of c. One can assume that it can be approximated by a white, 
delta-correlated noise term. This last term vanishes in correspondence of 
the absorbing states. 

If the mean field part /(c) converges towards an absorbing state, the 
dynamics is given by the competition between the diffusion and the noise 
term, and this produces the usual directed-percolation (or parity conserva- 
tion) phase transition, for which the presence of a stable, locally attracting 
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absorbing state is essential [9]. Far from the absorbing states (or if the 
absorbing states are not present, e > 0) the noise term is irrelevant. 

It is expected that a behavior more similar to that of mean-field can 
be observed if the diffusion term (or an equivalent mechanism) increases 
the spatial homogeneity. This is usually achieved in theoretical physics by 
increasing the dimensionality of the system, but this mechanism is unlikely 
to be observed in social networks. 

In real society, people are rarely arranged in a one-dimensional lattice. 
There are many proposed structures for social networks, but one feature, the 
small-world effect is generally present and its effect is that of an increasing 
spatial homogeneity. 

5. Smallworld bifurcations 

The Watts-Strogatz model [TU] is one of the simplest network models 
exhibiting the small-world effect which allows to smoothly change from a 
regular to a random lattice. We have therefore simulated the microscopic 
system on a regular one-dimensional lattice where a fraction p of links are 
rewired at random, and measured the behavior of the density c(t). Its 
asymptotic value, for large t, is denoted c(oo). 

In Fig. g we show the return map (c(t + 1) vs. c(i)) of the density of 
opinion 1 as results from the actual simulations, together with the mean- 
field predictions, for various values of p. One can see that for p = 0, the 
density simply fluctuates around its mean value, 0.5. By increasing the 
fraction of long-range links, the distribution of points approaches more and 
more the mean-field predictions. 

We show in Fig. [5] the bifurcation diagram obtained by increasing the 
probability of long-range connections, p, for different values of J. 

By comparing the plots of Fig. [2] and Fig. [5] one can see that one has 
close resemblance between the bifurcation induced by J and that induced by 
p. However, we were not able to find an analytic mapping. Notice that this 
hypothetical map is nonlinear: the first bifurcation occurs at larger values 
of p when decreasing J, but the second bifurcation moves in the opposite 
direction. 

For p above (roughly) 0.8, the distribution of points corresponds to that 
of mean-field, so one can assume that the small-world threshold for the 
chaotic oscillations of these models occurs at high values of the fraction of 
rewired links p, differently from the usual small- world effect. 

For J > the system is "ferromagnetic" (the transition probabilities 
are monotone with c, and so is c'(c). Therefore there are at most three 
fixed points corresponding to the intersection of the curve c'(c) with the 
bisectrix (one for c = 0.5 and two symmetrically placed). By increasing J 



p = 0.25 P = 0A 
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Fig. 4. Density of opinion 1 (c) from microscopic simulations (crosses) plotted 
as return map [c(t + 1) vs. c(t)] and mean-field predictions for various values of 
the long-range connection probability p. Here J = —3, k = 20, q = 0.1, e = 0.2, 
= 1000, transient of 1000 time steps and plot of 200 iterations 



a bifurcation appears (not reported here), very similar to that observed in 
an Ising model, but the branches corresponds to stable fixed points, not to 
cycles as for J < 0. This scenario is confirmed by numerical simulations. 
For p = 0, the system is ferromagnetic and finite-range, it is expected that 
this bifurcation is really stable only in the presence of absorbing states, i.e., 
for e = 0, while for other values of e it is only a metastable state (in real 
simulations). However, for p > the system is actually infinite-range and 
therefore a phase transition can be present also in one-dimensional systems 
with finite interaction terms. It might be that a transition from metastable 
to stable states occurs for some finite values of p, even though we suspect 
that the threshold is for p — > 0. We have not investigated this aspect in 
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Fig. 5. Small- world bifurcation diagram for J = — 2, . 
e = 0.2, N = 5 • 10 5 , transient of 2000 time steps and 2 different initial conditions 
plotted for 20 iterations. 



detail. 

6. Conclusions 

We have studied an opinion model that exhibits, at the mean-field level, 
a period-doubling cascade to chaotic oscillations, by varying the coupling 
parameter J. The observed quantity is the average density of opinion 1, c. 
Actual simulations on a one-dimensional lattice, in the absence of absorbing 
states or in the "active" phase, show microscopic chaos, i.e., incoherent local 
oscillations; so that the density c fluctuates around 0.5. 

By rewiring a fraction p of local connections to a random site, we trigger 
a small- world effect: the density c exhibits a bifurcation diagram that resem- 
bles that obtained by varying J. These small-world induced bifurcations are 
consistent with the general trend: long-range connections induce mean-field 
behavior. However, this is the first observation of such effect for a system ex- 
hibiting a chaotic mean-field behavior. Indeed, the small-world effect makes 
the system coherent (with varying degree). We think that this observation 
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may be useful since many theoretical studies of population behavior have 
been based on mean-field assumptions (differential equations), while actu- 
ally one should rather consider individuals, and therefore spatially-extended, 
microscopic simulations. The well-stirred assumption is often not sustain- 
able from the experimental point of view. However, it may well be that there 
is a small fraction of long-range interactions (or jumps), that might justify 
the small- wo rid effect. In particular, it would be extremely useful to derive 
a general "rescaling" formulation allowing the estimation of the effective 
value of parameters given a certain degree of small- world connections. 

For what concerns our specific opinion model, we can draw some socio- 
logical consequences from our simplified assumptions. We simulated social 
groups with "frustrated" behavior, i.e., conformistic for a strong social pres- 
sure and anticonformistic for "marginal" behavior like fashion or dressing. 
Such a scheme can be probably applied to many societies in the transition 
phase from traditional to non-traditional behavior, but also to social micro- 
cosmos in western societies, in which social norms are hardly broken. It is 
plausible that these frustrations may trigger oscillations, possibly chaotic. 
A social initiative promoting homogenization, or the social mixing due to 
living or working conditions, could act by favoring long-range interactions 
and triggering coherent oscillations. Such oscillations could be identified in 
the sudden explosion of violence or pathological trends (say suicide, self- 
mutilations, etc.). 
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